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Abstract. In spintronics, the ability to transport magnetic information often 
depends on the existence of a spin current traveling between two different magnetic 
objects acting as source and probe. A large fraction of this information never reaches 
the probe and is lost because the spin current tends to travel omni-directionally. 
We propose that a curved boundary between a gated and a non-gated region within 
graphene acts as an ideal lens for spin currents despite being entirely of non-magnetic 
nature. We show as a proof of concept that such lenses can be utilized to redirect 
the spin current that travels away from a source onto a focus region where a magnetic 
probe is located, saving a considerable fraction of the magnetic information that would 
be otherwise lost. 
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The field of spintronics is based on ingenious ways of controlling not only how 
the electron charge flows across nanoscaled materials, but also how the electron spin 
propagates within the same environment pp. One of the key challenges in the field is to 
find materials with large spin-coherence lengths, so that the information contained in 
the electron spin is not lost after a short propagation distance. Carbon-based materials 
such as nanotubes and graphene have been in the scientific limelight for over a decade 
now and among their many interesting properties is the fact that they have unusually 
long spin-coherence lengths (21 El HJ El E] . Therefore, it seems like a perfect match to 
study spintronics in such carbon-based materials, as confirmed by the growing number 
of studies on the subject [H El El HH HOI HH H21 H31 HS1 HB]- 

One of the peculiarities of graphene is the linear dispersion relation for electrons 
at the Fermi level region, which makes electronic wave packets travel dispersionlessly. 
This is contrary to what is commonly expected of typically dispersive quantum particles 
described by packets with wave-vector components that travel at different speeds. Such 
a lack of dispersion in graphene suggests that its electrons must evolve very much like 
typical classical waves. As a result, analogies with optical systems have been made in 
an attempt to widen the applicability of graphene [T71 HB] . Here we propose that this 
optical analogy be further explored to generate graphene-based lenses that are entirely 
non-magnetic and yet capable of concentrating spin currents onto a focal point. 

Spin currents correspond to a net flow of electron spins that may or may not 
coexist with charge current. One way of producing such currents is by pumping 
spins into a non- magnetic metal through the precession of an adjacent magnetization 
[19, 20J. In this case, angular momentum from the moving magnetization is transferred 
to the conduction electrons, creating a spin disturbance that propagates throughout the 
conducting material. It is noteworthy that spin currents may also be generated by a 
rotating magnetic field, as reported in reference [21J, with no need of magnetic atoms 
being attached to the graphene lattice structure. Let us emphasize that this spin flow is 
produced without an applied voltage and involves no net electrical current. Moreover, 
the traveling spin current can in principle be used to remotely excite other magnetic units 
located a long distance from the source of magnetic precession, so that the conducting 
material acts as a waveguide of magnetic information [221 [23]- The obvious drawback 
of using this method for transporting information is that because the initially localized 
precession tends to travel omni-directionally, a good fraction of this information is lost 
if it is not reabsorbed by another magnetic moment. Ideally, one would like to redirect 
the spin current and focus it towards a magnetic object capable of absorbing the spin- 
precession energy that would be otherwise lost. 

The ease with which graphene is gated suggests that it is straightforward to make 
its dispersion relation spatially dependent according to the gate shape. Because the 
gate strength defines the velocity with which electrons travel across a material, one 
can envisage a way in which spatially dependent refraction indices are engineered for 
electrons propagating in graphene. Both refraction and diffraction phenomena affect the 
propagation direction of a wave and are able to induce focussing patterns. In this work, 
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we show with a fully quantum mechanical approach how spin-current focussing may 
be achieved by dividing a large graphene region into two parts, one of which is under 
the action of a gate voltage. It is worth stressing the enormous benefit of controllably 
redirecting the spin-current in order to minimize the loss of magnetic information in 
spintronic devices. 

Let us consider a graphene sheet with a few magnetic atoms substitutionally 
inserted into its hexagonal atomic lattice. Magnetic atoms, usually transition metal 
elements, are known to display sizable magnetic moments [2*1] 1251 126] when inserted 
in graphene and will be used here as the source of spin current. Other objects such 
as magnetic nanoparticles, ferromagnetic substrates or even vacancies [27] can also be 
used but substitutional atoms are by far the simplest and yet capture the essence of the 
phenomenon we intend to describe. We also consider one stripe of magnetic atoms as 
schematically illustrated in figure [Tj While a pristine stripe of substitutional impurities 
is unlikely to exist, here it serves the purpose of (a) testing how the system responds 
to the size of the magnetic material responsible for generating the spin current and 
(b) inducing spin currents that quickly develop into plane waves. Furthermore, we 
assume the system to have a curved boundary separating two different media, one of 
which under the action of an external gate voltage V g which effectively controls the 
electron density in the region. Figure 1(a) shows a schematic diagram of the system 
under consideration. Note that the boundary between the gated (gray) and non-gated 
(white) regions is curved and chosen to be circular with a radius R. The obvious 
problem of considering such a setup is the difficulty of describing this system without 
paying incredibly hefty computational prices given the enormous number of atomic sites 
required. For that reason we choose to work with finite-width sections of graphene with 
lateral periodic boundary conditions (PBC), where the repeated unit cell is the part of 
figure pja) delimited by the two horizontal dotted lines. Note that flat boundaries are 
also straightforwardly considered by taking the limit R — > oo, as shown in figure Pub). 




As previously mentioned, the pumped spin current is generated by the precession 
of a magnetization in contact with a conducting medium. We assume that the 
magnetization is originally in equilibrium, pointing along an arbitrary z-direction, and 
that it is set into precession by a time-dependent transverse field h±(t). In practical 
terms, there are different ways of producing as well as probing this time-dependent 
perturbation [28j 12S1 EQ] . To determine how graphene transports spin currents, we must 
assess how such a localized magnetic excitation propagates across the structure. When 
such a spin current flows across the normal metal it creates a time dependent local 
spin disturbance that may be described by the transverse spin susceptibility. If the 
amplitude of h±(t) is sufficiently small, one may use linear response theory to relate the 
spin disturbance at site m with the driving field applied at site £. The relation for a 
harmonic perturbation with amplitude h is given by 5(S+(i)) = gfiBh e~ luJt Xm~e{ UJ )/'2, 
where = S^+iS^, and S^ y represent the local transverse spin components at site m. 
It follows that \5{S^{t))\ = gfj, B h \x^}{(^)\/^ and, therefore, the off-diagonal transverse 
spin susceptibility provides the amplitude of the spin disturbance caused by the spin 
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Figure 1. Schematic diagrams representing the structures under study, (a) An 
infinitely large graphene sheet is under the effect of a gate voltage that acts only 
in a limited region of space (represented by the gray area). The gated region has a 
boundary with a curvature defined by the radius R that lies at a distance L\ from the 
the magnetic atoms. The section delimited by the two horizontal dotted lines is the 
unit cell used in representing the system with lateral periodic boundary conditions, (b) 
Flat boundaries are easily accounted for in the limit R — > oo. Dashed lines represent 
the locations where the spin disturbance will be probed. 

current that flows through the normal metal. Here, we investigate how this excitation 
disturbs the spin balance of the system not only where the magnetic atoms are located 
but also, and more importantly, how the local spin dynamics is affected within graphene. 

To calculate the spin susceptibility j2D E2] one needs the Hamiltonian describing 
the electronic structure of the unperturbed system, which we assume is given by 
H = claCja + He, a ^.^la^la + H z . Here, jij represents the electron hopping 

between nearest neighbor sites i and j, c ifJ creates an electron with spin a at site i, the 
sum in £ is over the sites occupied by magnetic atoms, n i(T = & ta c.fo is the corresponding 
electronic occupation number operator, and U represents an effective on-site interaction 
between electrons in the magnetic sites, which is neglected elsewhere. With the lateral 
PBC imposed, the Hamiltonian matrix is considerably reduced. Finally, Hz plays the 
role of a local Zeeman interaction that defines the z-axis as the equilibrium direction 
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Figure 2. Magnitude of the spin susceptibility calculated at resonance frequency for 
different positions along the horizontal dashed lines shown in figurc[I] The insets depict 
the spin susceptibility probed at the carbon sites situated along a vertical line located 
at position 5 of the horizontal axis shown in figure [I] The dotted lines in figures (a) 
and (b) mark the position of these vertical lines. All graphs are for V g = 2.4 eV. (a) 
Results for a flat boundary (R — > oo). (b) Results for a curved boundary (R = 3a) 

of the magnetization. The Hamiltonian parameters can be obtained from density- 
functional-theory calculations so that the electronic structure of the doped system is well 
described [2SIES]- Although the presented results are for Mn atoms, other substitutional 
magnetic impurities may be employed. In our calculations we fix the Fermi energy 
Ep = 0, and use jij = 2.7 eV for the nearest-neighbour hopping terms. We take the 
number of <i-electrons in the Mn sites n m = 1, U = 20 eV, and assume a local Zeeman 
energy splitting of 1 meV, which is of the order of magnitude of magnetic anisotropy 
energies found for transition metal atoms adsorbed to metallic surfaces |35J GUI- For 
simplicity we consider the electronic hopping between carbon atoms to be the same as 
the hopping between carbon and the magnetic impurity. Within a five-fold degenerate 
orbital scheme, this requires a large value of U to obtain a reasonable value for the Mn 
magnetic moment. Spin-orbit coupling, which has been shown to play a role in the case 
of magnetic adsorbates[40j, is neglected in the substitutional impurity case for being 
very small compared to the other relevant energy scales. 

The time-dependent transverse spin susceptibility is defined as Xm/if) = 
— lQ(t)([S+(t), Sl(0)]), where 6(x) is the heaviside step function, and and S m 
are the spin raising and lowering operators at site m, respectively. The indices i 
and m refer to the locations where the field is applied and where the response is 
measured, respectively. In our case, we induce a precession of the magnetic moments 
on sites £ and we wish to probe the resulting spin disturbance at an arbitrary site m. 
Note that the probing site may or may not contain another magnetic object without 
affecting the conclusions here obtained. The effect of a magnetic object occupying 
the probing position simply amplifies the resulting spin disturbance. In what follows 



Graphene as a non-magnetic spin-current lens 



6 



we present results for which there are no magnetic sites on the probing region. In 
this case, this precession response is fully described by J2eXm/(t)- Within the random 
phase approximation, this susceptibility may be expressed in frequency domain, and in 
matrix form it is given by = [1 + U]" 1 x (a>), where x° is the Hartree-Fock 

susceptibility given by 

xl» = 7T- [ du'f{u') {[g\ m {uj') - g^{u')}gi e {u' + u) 

271 J-oo k 

+[gU"') - 9 m \{u'Ml{u' - u)} . (l) 

Here, g^ei. 00 ) an d 9me( u ) represent the time Fourier transforms of the retarded and 
advanced single-particle propagators for an electron with spin o between sites m and £, 
respectively, and f(u) is the Fermi-Dirac distribution function. Following reference |32j, 
it is possible to extend some those integrals to the complex plane, thereby avoiding a 
great deal of numerical difficulties. The one electron propagators are obtained by a well- 
established numerical methods IM]. x° describes single-particle (Stoner) excitations 
whereas the enhanced susceptibility x, calculated within the RPA, takes into account 
correlated electron- hole pairs, that correspond, in an extended ferromagnet, to collective 
modes known as spin waves [35] . 

We have recently studied the spin disturbance in carbon nanotubes as a function of 
time and have identified that pulsed magnetic excitations travel across these materials 
with their corresponding Fermi velocities, with very little deformation and with 
adjustable levels of attenuation, something that has been attributed to the distinctive 
linear dispersion relation displayed by nanotubes [22]. Graphene is no different in this 
regard and, for having similar linearities in its dispersion relation, should behave in 
exactly the same way. In fact, we have demonstrated this similarity with graphene 
ribbons by showing that localized spin excitations also propagate dispersionlessly and 
with negligible loss. Furthermore, we have also argued that a gated region between the 
magnetic object generating the spin current and the locations where the spin disturbance 
is probed may be capable of modulating the spin precession between on- and off-states, 
giving rise to a transistor behavior for the pumped spin current that functions without 
the need of a bias to drive the current [23J. 

Rather than studying the time-resolved spin excitations [22], in this paper we focus 
on the frequency- dependent response. The key quantity for us is therefore |x m ^(w)| 
because it is proportional to the amplitude of the spin disturbance at site m due to a 
time-dependent transverse magnetic field applied at site I. In other words, it corresponds 
to the precession amplitude acquired by the electron spin localized at an arbitrary site 
m due to a harmonic perturbation of frequency u produced at the magnetic site t. 
Unsurprisingly, when investigated as a function of frequency this quantity displays a 
very distinctive peak close to its Larmor-frequency resonance. Without any loss of 
generality, we use this value as our choice of excitation frequency. 

Figure 2(a) shows the results for |x m ^(w)| probed at different sites m, all of which 
following the horizontal dashed line depicted in Figure 1(b). This result corresponds to a 
flat interface (R = oo) separating the two different media, one of which under the action 
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Figure 3. Color-coded plot of the spin susceptibility as a function of the position 
within the gated region depicted by the gray areas in figure [l] The Top (Bottom) panel 
corresponds to a flat (curved) boundary with parameters specified at the captions of 
figure [2j 



of a gate voltage V g = 2.4 eV. There is a noticeable oscillation in the amplitude of the 
spin precession as a function of the probing position along the axis of the ribbon, which 
can be ascribed to the inevitably oscillatory spin polarization induced by the presence of 
magnetic objects in contact to conduction electrons. As previously mentioned, because 
the precession amplitude is comparatively increased for larger spin moments the same 
oscillation appears in the susceptibility results. The inset of figure |2](a) shows that 
the response is fairly homogeneous when probed along a direction perpendicular to the 
ribbon axis, located at position 5 of the horizontal axis shown in figure [T|b). This means 
that the spin precession caused by the spin current travels into the gated medium as a 
plane wave without favoring any subpart of the system. 

A different response in the precession amplitude is found in figure [2^b) when the 
flat interface is replaced with a curved surface of radius R = 3a, where a is the lattice 
parameter of graphene[41j. The gate voltage is the same as the one used in figure 
[2^a) but the amplitude of the oscillations is nearly twice as large. The oscillatory spin 
polarization induced by the proximity to a magnetic object is known to be primarily 
dependent on the details of the medium in which the polarization occurs, indicating that 
a change in the interface shape cannot enhance the induced moments very significantly. 
Therefore, the increase in the amplitude of |x m ,^(w)| is likely to result from a focussing 
process. A more convincing evidence that the spin current is being focussed by the 
curved interface is shown in the inset of figure [2^b), which displays |x m /(o;)| probed on 
a line that is perpendicular to the ribbon axis, located at position 5 of the horizontal 
axis shown in figure [l|a). A maximum is clearly seen at the central region of the 
ribbon, indicating that the peaks seen as a function of the axial positions are not 
distributed homogeneously across the ribbon width but are concentrated at a small 
region of space acting as a focus for the spin current. This is even more conclusively 
shown in figure [3j where 2-dimensonal color plots representing the precession amplitude 
as a function of the probing position in the graphene ribbon is displayed. The upper 
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panel corresponds to the flat interface and is rather featureless. That confirms our 
earlier statement that the precession amplitude is rather homogeneous for flat surfaces. 
For the case of curved interfaces, however, the bottom panel of figure [3] shows a region 
where the spin susceptibility reaches a localized maximum acting as a center where the 
spin current is focussed. It is worth noting that this region is surrounded by other areas 
where the spin susceptibility is very low. This is understood when we compare the 
bottom and upper panels of figure [3j indicating that the precession energy that would 
be homogeneously spread across the ribbon width in the flat-surface case is redirected 
towards a focus when in the presence of a curved interface. As a result, a large part of 
the magnetic information contained in the precession motion of the electron spins would 
not be probed and would be completely lost. By playing with the combination of gate 
voltage and geometry of the interface separating two parts of a graphene sheet, we can 
minimize this information loss. Finally, it is worth stressing that the results reported 
here are very robust and remain qualitatively the same for a wide range of diameters 
and voltages. This opens a number of possibilities for using the energy contained in 
the precession of magnetic moment as a source of magnetic information. The fact that 
lenses can be used to concentrate the spin current onto a small region of space is ideal 
to amplify spin signals that could be otherwise very weak. 

In summary, we propose that the spin current originated from a precessing 
magnetization in contact to a graphene sheet can be controllably redirected towards 
a focus region, which could reduce the loss of magnetic information stored in the spin 
dynamics of such systems. Furthermore, by probing the spin current at the focus region 
one is likely to enhance the sensitivity of any spin pumping device of this type. 
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